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Abstract 

The effect of anisotropy in the nearest-neighbor spin interactions that couple N > 2 
consecutive spin-1/2 antiferromagnetic chains is studied theoreticaUy by considering the 
limit where the coupling is purely of the Ising type. An analysis based on the equivalent 
Luttinger model reveals that the groundstate is an Ising antiferromagnet in general. 
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The recent discovery of the "ladder" materials spawned from research in high-temperature| 
superconductivity has renewed interest in the physics of coupled antiferromagnetic chains.-^ 
The former systems are composed of magnetically isolated ladders of spin- 1/2 moments 
that experience an effective exchange interaction between all nearest neighbors. Both ex- 
periment and theory find that the S = 1/2 antiferromagnetic ladder shows a spin gap 
of order the exchange coupling constant between chains. ^'^ This result, however, should 
be compared with that corresponding to a single chain, which exhibits no spin gap.^ In 
general, an analysis of the weak and of the strong-coupling limits reveals that no spin gap 
appears for an odd number of chains. 

Deviations from isotropy in the Heisenberg spin coupling that results from the ex- 
change interaction can occur naturally, however, and have interesting theoretical conse- 
quences. It is known, for example, that Ising anisotropy produces a spin gap in the case of 
a single chain.^"^ In the present context of coupled antiferromagnetic chains, the following 
question then arises: What effect does anisotropy in the Heisenberg spin couplings that 
run perpendicular to the spin-1/2 chains have on the ground state? In response to this 
query, we shall study here the low-temperature physics of a finite number N of consecutive 
spin-1/2 XX Z chains^ coupled via a nearest-neighbor Ising interaction. The calculational 
strategy will be to first transcribe the problem to that of interacting spinless fermions 
utilizing the Jordan- Wigner transformation, and to then apply the well-known abelian 
bosonization technique. ^'^ On this basis, we arrive at the following conclusions valid for 
N > 2 Ising-coupled antiferromagnetic chains: (i) The groundstate is in a pinned spin- 
density wave (SDW) state that is commensurate with the lattice,^ and that exhibits a spin 
gap; (ii) The SDW slides along the direction parallel to the chains in the case where it is 
incommensurate,^ which can result from the application of an external magnetic field. It 
is noteworthy that the bosonization analysis employed here is closely related to a fermion 
analogy for the Lawrence-Doniach model description of layered superconductors. ^-"^'^^ We 
now pass to the calculations. 

The Hamiltonian for N Ising-coupled antiferromagnetic chains that are sequentially 
ordered can be divided into parallel and perpendicular parts, H = + H±, where 

AT 

1=1 i 
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and 

N-l 

-f^-L = X] Jz^liSli+i (2) 
1=1 i 

describe respectively the spin couplings within and in between the corresponding XXZ 
chains. Here, the spin operator Si^i acts on spin-1/2 states lying at the i^^ site of chain I, 
and hz denotes an external magnetic field directed along the z direction. We presume an 
antiferromagnetic sign, j|y > 0, for the intra-chain XY coupling, and we shall set h = 1 
throughout. The Jordan- Wigner transformation^'^ then yields a system of ideal spinless 
fermions, with a degenerate energy spectrum = — jj^cos/ca as a function of momentum 
k along each chain I, that interact through both the intra-chain and the inter-chain Ising 
couplings, j| and J^. Here, a denotes the parallel lattice constant. In this language, 
an antiferromagnetic chain then corresponds to a half-filled band, Sk < 0. Note that 
inter-chain XY coupling introduces unwieldy "string" contributions^ into the Hamiltonian, 
which is one reason why it has been omitted. The continuum limit can be taken d la Kogut 
and Susskind under these conditions,^ and we thereby obtain the following Luttinger model 
for the parallel and perpendicular pieces of the Ising-coupled antiferromagnetic chains: 



Af „ r 

H\\=J2 4ya(^L}idM - Rjid^Ri^ + UlaL]R]RiLi + h,{L\Li + r]ri) 



1=1 

and iy_L = iy_L,i + -f^_L,3 + -f^_L,4, with a backscattering term 

N-l 



(3) 



1=1 

an inter-chain umklapp term 

N-l 



H^,i = J2 f ^^-^^^^ [LjRj+^Li+iRi + H.c] , (4) 

/ -I ^ 
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1=1 

and an inter-chain forward scattering term 

N-l 



H±,4 = J2 [ ^^"^^^^ • + ■■■■ (4+1^1+1 + Rj+iRi+i) ■■ ■ (6) 
1=1 

Here e~'^'^''^ Ri{x) and e'^'"^^Li{x) denote field operators for right and left moving spinless 
fermions that move along the l^^ chain, with a Fermi surface at ifci?, while the symbols 



': :' represent normal ordering. ^'^ The intra-chain umklapp term^'^ that has been omitted 
from Eq. (3) will be discussed at the end of the paper [see Eq. (15)]. 

We first consider two chains, in which case the above Luttinger model can be treated 
exactly. This is achieved by observing that the chain index can be interpreted as a pseudo 
spin label. Equations (3)- (6) then describe the Luther-Emery model for pseudo spin- 
1/2 fermions in this instance. Since such fermions experience pseudo spin-charge 

separation, we have that the coupled chains factorize following i^n -|- H± = Hp + iy^-, 
where 

Hp=27rVpJ2 Yl pM)pA-(1) + 9pYpR{(l)pL{-q) + H^^^ (7) 

q>Qj=R,L q 

=27TVa Y +9ctY ^Ri(l)^L{-q) + -H"±,i (8) 

q>Oj=R,L q 

are the respective commuting portions of the Hamiltonian. Here, Pj(q) = 2~^/^[pj(q, 1) -|- 
Pj(q,2)] and aj(q) = 2~^/'^[pj(q,2) — pj(q,l)] are the particle-hole operators for pseudo- 
charge and pseudo-spin excitations, with pR(q,l) = X^fc '^J + ^)'^'(^) Pl{q,1) = 
J2k^}(^'^ /c)6i(/c). The operators ai{k) and 6/(/c) respectively annihilate right and left 
moving electrons of momentum k on the l^^ chain. Also, the Fermi velocities and interac- 
tion strengths for each component are renormalized by the inter-chain forward scattering 
process (6) to 

vp,^ -Jj, ± J^/27:, (9) 
gp,,=U^}±J^, (10) 

where the +(— ) signs above correspond to the p{a) label. Application of the bosonic 
representation^ 

Ri{x) ^(27ra)-i/2exp[z(47r)i/2(/,fl(a;, 0] (11) 
Li{x) ^{2na)-^/^exp[-i{4ny/^(t>L{x, I)] (12) 

for the spinless fermions, where (t)j{x^l) — limQ,^o27rL~^ g~''^exp(— — iqx)pj{qj) 
are the bosonic fields corresponding to right and left moving fermions {j = R,L), reveals 
that the spectrum of the pseudo-spin sector (8) has a gap, A^. ^ 0, for —g^ < |J^|.^^ 
(In general, ~ is the momentum cutoff of the Luttinger model.) At half-filling 
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{hg = 0), where Akp = 27r/a and the umklapp process (5) is at work, similar considerations 
indicate that a gap, Ap ^ 0, opens in the pseudo-charge spectrum (7) for —gp < \J^\. 
Comparison of Eq. (10) with these conditions yields the phase diagram shown in Fig. 1 for 
two Ising-coupled spin-1/2 XX Z chains. We remind the reader that the spinless fermions 
corresponding to the pseudo-charge (7) and the pseudo-spin (8) systems are noninteracting 
along the respective Luther-Emery lines Qp = GirVp/b and = 67rf(j/5, at which point the 
gaps have value Ap^^- = the other hand, the Coulomb gas analogy-*^^ as 

well as mean-field theory^^ indicate that the latter vanish exponentially as the Luttinger- 
liquid state is approached; e.g., Ap^. ~ | J^|exp[— 27rVp^o./(|<^i"| + 9p,(t)]- 

We now address the issue of the physical character of the two-chain system in the 
case where all spin coupling are antiferromagnetic, which means that Ap 7^ (see Fig. 
1). Since the longitudinal paramagnetic susceptibility is simply given by the pseudo- 
charge compressibility, we have that this quantity follows the activated behavior Xz oc 
exp(— 2Ap//cbT) at low temperatures. Also, pseudo charge-spin separation implies that 
the specific heat is given by the sum Cy = Cp + Co- of the respective charge and spin 
contributions, each of which follow the activated behavior Cp^a exp(— 2Ap^o-/A;bT) at 
low temperatures. The question of long-range order at zero temperature can be attacked 
with the bosonization method (11) and (12). Following Luther and Peschel,^ the static 
transverse spin correlator on the same chain may be calculated using the formula 

that is valid to lowest order in e*^^^^. Pseudo spin-charge separation then implies that 
this correlator has the asymptotic form ('S'^i'S'^^i) = e^'^^^^G^r^g{x)Gj^g{x)/ Aira, where 
G^g{x) = {a/x)^''e~^/^'' is the auto-correlation function for pseudo-triplet superconductivity,^ 
with Cooper pairs of effective unit charge, and where Gifg{x) is equal to the previous mod- 
ulo the symbolic replacement p a. Here, 9a,p = 1/4 and a/i(T,p = Ao-,p/fo-,p- Hence, we 
arrive at the result {S^^-^Sq ^) ~ e*^'^^^(a/x)^/^e~^/^'"f , where ^xy/a = {^p/vp + /v„)~^ 
is the (finite) XY correlation length that signals short-range transverse spin correlations. 
Lorentz invariance in the dynamical correlator G'>f^\x,t) for pseudo-triplet superconduc- 
tivity then implies that the dynamical transverse spin correlator (13) exhibits a spin gap 
Axy = Ap + Ao-. Similar calculations reveal that the longitudinal spin correlator has the 
asymptotic form (S^ -^^S^ -i) = cos{2kFx)G''^ly^,{x)G^gl,y^,{x)/{4:7^a)^, where G^gl,y^,{x) 1 



is the autocorrelator for pseudo SDW order/ and where GsDwi^) obtained again 
through a trivial symbohc replacement. We therefore find that such Ising-coupled an- 
tiferromagnetic chains display strict long-range order of the Ising type, like that displayed 
by an isolated Heisenberg chain in the presence of Ising anisotropy.^"^ (See Table I for 
a listing of the relevant static correlation exponents). In conclusion, the system is in a 
pinned SDW state that is commensurate with the lattice, and that necessarily exhibits a 
spin gap. 

The above results imply that this SDW state is incommensurate^ [kp ^ 7r/2a) for 
large enough Zeeman energy splittings^*^ \hz\ 3> Ap, in which case umklapp processes 
become irrelevant. Then by Eq. (7), the pseudo-charge sector is in a pure Luttinger 
liquid state, Ap — > 0. Also, since the bosonic field 0p = + ^2) now represents a 

Goldstone mode, where (j)i{x) = (j)L{xi I) + (f>R{x, I), we have that the SDW can freely slide 
along the chain direction. The list of physical properties that we outlined above must then 
be revised as follows. Since the pseudo-charge system is now compressible, the longitudinal 
paramagnetic susceptibility no longer shows a spin gap. In particular, we have that Xz = 
i^Jxy + Jz + ^J^)~^ a-t zero temperature. Likewise, the low-temperature specific heat is 
now dominated by the pseudo-charge component Cp oc T. Last, the sliding Goldstone mode 
results in only algebraic longitudinal spin correlations {S^^iSq^i) ~ cos{2kFx){a/x)^'' , 
with exponent Kp = {27rvp — QpY^'^ / {2'KVp + QpY^^- Note that here refers to the 
deviation of the magnetization with respect to it's average value, Xzhz- And although 
the transverse spin correlations remain finite in the the present incommensurate case, the 
XY spin-correlation length, ^xy/*^ — ^al is now larger. 

We shall now treat the general case of N > 2 chains. The previous results, (9) and 
(10), obtained for the case of two Ising-coupled chains indicate that the inter-chain forward 
scattering process (6) can be neglected in the limit of weak coupling, \J^\ <^ jjy, | 
This shall be assumed throughout in the present discussion. To begin with, we shall also 
neglect all umklapp processes. The Luttinger model for the Ising-coupled antiferromagnetic 
chains in such case reduces to the sum of Eqs. (3) and (4), which describes a generalized 
backscattering model. A mean-field analysis^^ of this model for large finds that long- 
range order of the charge-density wave type, {Lj{x)Ri{x)) oc exp[z(47r)^/^0/(a;)], is stable 
for \J^\ > — 2j|. Notice that this agrees with the exact results for two chains (see Fig. 1). 
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In addition, the application of the abehan bosonization technique, (11) and (12), yields 
the action 

r N N-i >| 

5ld = W dxodx^ I Yl - ^o"' E cos[(47rK) V2(0;^^ _ 0^)] K (m) 

U=l 1=1 J 

related to the Lawrence-Doniach model of layered superconductivity,^^ where (j)'i{x,t) rep- 
resents the time evolution of the bosonic field operator, (f)i{x), after renormalization.^ Here, 
K = e^'^ is one-half the exponent for Ising order, with the angle i/j set by the condition 
tanh2'i/' = — 2j|/7rjjy, while = {iv'pt^x) is a two- vector, with a renormalized Fermi 
velocity v'p — j|j,a scch2'i/;. Also, oc (J^)~^/^ denotes the XY spin-correlation length. 
An analysis based on the equivalence of (14) to a layered Coulomb gas^^ indicates that 
this term is relevant in the renormalization group sense for \ J^\> — 4j|. Notice that this 
condition is roughly consistent with the phase diagram corresponding to two chains (see 
Fig. 1), as well as with the mean-field result just cited. Also, it is clear from (14) that 
collective sliding along the chain direction, (t)i{x) (j)i{x) + Scpp, represents a Goldstone 
mode. All of these pieces of information when put together then strongly suggest that 
the groundstate of weakly Ising-coupled chains corresponds to an unpinned SDW state if 
umklapp processes are absent. 

Suppose now that we include the umklapp processes present (at half-filling) in zero 
field. The addition of the inter-chain process (5) as well as the intra-chain one"^"^ then 
leads to the bosonic action 



N-l 



^ Y 2^o-'cos[(47rK)V20|^^]cos[(47rK)V2<^|] , (15) 
1=1 ) 

where ^2 represents the XY spin-correlation length of an isolated XXZ chain. If we then 

fix a particular chain I and integrate out the neighboring fields, -"^^ we obtain the effective 

action 

Si = i^ d^x^-{d^4>'i? - er'cos[(47rK)VVj] - e2-'cos[2(47rK)VVj]}, (16) 

which is a sine-Gordon model^^ with a first harmonic due to the intra-chain umklapp pro- 
cess. Note that the coefficient ^ of the base term is not free, it being proportional to 
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the appropriate average (cos[(47ri^)-^/^0^_j^]) + (cos[(47ri^)-^/^0^^.;^]) that results from the 
integration. Within the equivalent Coulomb gas description, this inter-chain term corre- 
sponds to unit charges, while the former intra-chain term proportional to $,2'^ corresponds 
to double charges. Since double charges in general dissociate into unit charges if the system 
is in the (spin-gap) plasma phase, we conclude that the first harmonic term is irrelevant 
in the presence of the base term. In the limit of weak coupling, — > 0, this means 
that a spin-gap (^1 < 00) opens for j| > — (Stt/IO) jjy, and that bound states appear 
for j| > O.-*^^ In contrast, a spin gap will open only in the presence of Ising anisotropy, 
Jz > Jxy, for the case of an isolated XX Z chain (J^ = 0).^ On the basis of what is known 
for a single XX Z chain,^'^ we conclude that Ising coupled spin-1/2 XX Z chains with 
all-antiferromagnetic couplings are generally in an Ising-antiferromagnetic state. Notice 
that this claim relies on the validity of the reduction (16) to a single spin chain, which 
clearly fails in the case of an odd number of Heisenberg-coupled chains. Nevertheless, 
the fact that we obtain a sliding SDW state in the absence of umklapp processes, coupled 
with the fact that umklapp processes only reinforce (commensurate) SDW order, strongly 
indicates that such a reduction is indeed correct. 

To conclude, we find that anisotropy in the spin-coupling that may exist in between 
antiferromagnetic chains can dramatically change the nature of the groundstate of each 
chain. In particular, the present Luttinger model based analysis demonstrates that consec- 
utive spin-1/2 antiferromagnetic chains flow to the Ising antiferromagnetic flxed point in 
the presence of any amount of Ising coupling in between nearest-neighbor chains (see Fig. 
1). This contrasts with the case of a single chain, where the Ising antiferromagnetic state 
exist only in the presence of Ising anisotropy, Jz > Jxy- It also contrasts with the situation 
where all couplings are Heisenberg, in which case no spin gap exists for an odd number 
of chains. The latter suggests that a transition between the XY and the Ising antiferro- 
magnetic fixed points occurs at some intermediate perpendicular anisotropy, J^y ~ J^, in 
such case. 
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Figure Caption 



1. Shown is the phase diagram for two spin-1/2 XX Z chains coupled via a weak Ising 
interaction {J^)- The Ising antiferromagnetic (AF) regions {^p,a 7^ 0) are charac- 
terized by strict long-range spin order along each chain, while the XY AF regions 
(Ap,^ = 0) show dominant XY spin correlations (S^^iSq i) oc e*2'=J-*(a/a;)(^p '+'^- '^/^ 
along each chain I [see Eq. (13)]. The sign of the intra-chain Ising coupling (j|) 
switches between these two phases by passing through an intermediate XY dimer 
phase {Ap 7^ 0, A^r = 0) characterized by dominant inter-chain dimer correlations 
{S^iS~^2^Q 2^Q i) oc (a/x)^'^ for antiferromagnetic coupling > 0. Quadrapolar 
XY spin correlations {S^iS^2'So2'^oi) {oi/x)^p dominate in the opposing region 
(Ap = 0, Ap. ^ 0) at ferromagnetic couplings. Note that the superscript indicates 
local ferromagnetic order in between chains. 
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TABLE I. Listed is the correlation exponent r] obtained via the bosonization technique^ 
for various order parameters, 0{x), in the spin ladder; i.e., (O(x)Ot(O)) oc (a/x)"^. An- 
tiferromagnetic Ising coupling in between chains is assumed (see Fig. 1). Note that the 
value r/ = indicates strict long-range order, while rj = oo indicates short-range order. 
Below, we have Kp^^ = {2T^Vp^a - gp,aY/'^ / {^T^Vp^a + gp,aY''^ ■ 



Order Parameter 


T] {XY AF) 


ri {XY Dimer) 


rj (Ising AF) 


qv 

cz 
c+ c- 


Kp + K^ 


oo 
oo 


oo 


oo 
oo 



12 



(LADDER) 



\ XY 

\dimer 

AT Ah \ 


1 1 AC 

IbINo Ar 




XY AF<^' / 

/xY 
/ QUAD. 


Jz" 

ISING AF<'' 



